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Abstract. A number of observed phenomena associated with individual neutron star systems or neutron star popu¬ 
lations find explanations in models in which the neutron star crust plays an important role. We review recent work 
examining the sensitivity to the slope of the symmetry energy L of such models, and constraints extracted on L from 
confronting them with observations. We focus on six sets of observations and proposed explanations: (i) The cooling 
rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear 
pasta regions of the inner crust, (ii) the upper limit of the observed periods of young X-ray pulsars, confronting mod¬ 
els of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer, (iii) glitches from 
the Vela pulsar, confronting the paradigm that they arise due to a sudden re-coupling of the crustal neutron superfiuid 
to the crustal lattice after a period during which they were decoupled due to vortex pinning, (iv) The frequencies of 
quasi-periodic oscillations in the X-ray tail of light curves from giant fiares from soft gamma-ray repeaters, confronting 
models of torsional crust oscillations, (v) the upper limit on the frequency to which millisecond pulsars can be spun-up 
due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold fre¬ 
quency determined in part by the viscous dissipation timescale at the crust-core boundary, and (vi) the observations of 
precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering 
caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star 
just before merger. 

PACS. 97.60.Jd Neutron stars - 26.60.Gj Neutron stars, crust - 26.60.Kp Neutron stars, equations of state - 26.60.-c 
Neutron stars, nuclear matter aspects of - 97.60.Gb Pulsars - 21.65.Ef Symmetry energy 


1 Introduction 

The problem of constraining the isospin-dependence of nuclear 
forces is a vigorous field bringing together nuclear physicists 
and astrophysicists in theoretical, experimental and observa¬ 
tional studies. The motivations involve obtaining a better un¬ 
derstanding of the structural and dynamical properties of neu¬ 
tron rich nuclei and neutron stars, and an improved microscopic 
understanding of the nucleon-nucleon interaction and its emer¬ 
gence as a residual strong interaction of QCD. 

The interior structure of a neutron star (NS) is poorly con¬ 
strained because of uncertainties in the high-density behavior 
of the nuclear equation-of-state (EOS) and the plethora of ex¬ 
otic phases that have been theoretically postulated to exist to¬ 
wards the center of the star. Direct observational probes of the 
core of the neutron star are limited to neutrino and gravitational 
wave emission, both of which are yet to be directly detected. 
The outer regions of the star are more readily accessible via 
electromagnetic radiation from the surface and magnetosphere, 
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and provide probes of the physics of the neutron star crust and 
the EOS around nuclear saturation density. Such constraints 
are intrinsically interesting and help scaffold our understand¬ 
ing of the deeper layers of the star. In this review, we focus 
on a number of recently explored neutron star observables that 
probe crustal properties, focussing on their potential to provide 
constraints for the nuclear symmetry energy, the most uncertain 
part of the nuclear EOS at densities close to nuclear saturation. 

The isospin dependence of nuclear forces is manifest in 
dense nuclear matter as an energy cost of decreasing the proton 
fraction away from symmetric matter N = Z, a. cost that is en¬ 
coded in the symmetry energy as a function of density. There 
are many different notations for the symmetry energy and its 
expansion parameters around saturation density; let us intro¬ 
duce the one we use now. 

The degree of isospin asymmetry in a system may be ex¬ 
pressed locally in terms of the proton rip and neutron rin num¬ 
ber densities (which might vary in space) via the local asym¬ 
metry parameter 6 = — np)/n where n = nn + np is 

the total baryon number density, or globally via the parameter 
/ = {N—Z)IA, where W, Z and A are the total number of neu- 
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Fig. 1. (Color online) Recent constraints on the slope of the symmetry energy L from analysis of terrestrial experiments and astrophysical 
observations, some of which can be found summarized in this issue. Taken from ITol . The community average L ^ 60 MeV should be taken 
only as guide to the favored values of L that emerge from the wide variety of experimental evidence. 


trons, protons and nucleons in the system. For uniform nuclear 
matter, both parameters 6 and / are identical. Nuclear matter 
with equal numbers of neutrons and protons (6 = 0) is referred 
to as symmetric nuclear matter (SNM); nuclear matter with 
S = 1 is naturally referred to as pure neutron matter (PNM). 
Nuclei on Earth contain closely symmetric nuclear matter at 
densities close to nuclear saturation density po ^ 2.6 x 10^^ 
g cm“^ = 0.16 fm“^ = no. Nuclear experiments tend to con¬ 
strain the behavior of the binding energy of symmetric nuclear 
matter, E{n ^ nQ,6 ^ 0) to within relatively tight ranges, 
but direct ab initio calculations there are extremely difficult. In 
contrast, experimental data directly probing E{n no, 6 1) 

is impossible, but state-of-the-art ab initio calculations of PNM 
have led to significant constraints on E{n < no, ^ 1). By ex¬ 
panding E{n,5) about ^ = 0, 

E{n,5) = Eo{n) + S{n)5'^ + ..., (1) 


we can define a useful quantity called the symmetry energy 


S{n) 


ld‘^E{n,5) 

2 


( 2 ) 


which encodes the change in the energy per particle of nuclear 
matter as one moves away from isospin symmetry. This allows 
extrapolation to the highly isospin asymmetric conditions in 
neutron stars. The simplest such extrapolation, referred to as 
the parabolic approximation (PA) gives the relation 

^pnm(^) ~ Eo{n) E S{n). (3) 

Since our experimental constraints are dominated by results 
from densities close to no, it is customary to expand the sym¬ 
metry energy about X = 0 where x = obtaining 


S{n) — J E Lx + ^FCgyinX^ + •••, (4) 


where J, L and i^sym are the symmetry energy, its slope and 
its curvature at saturation density. The true values of the higher 
order symmetry energy parameters L, ATsym, ••• are still some¬ 
what uncertain, and the measurement of L in particular has 
been the subject of an intense experimental campaign by the 
nuclear physics community in recent years using nuclear probes 
such as masses, neutron skins, nuclear electric dipole polariz¬ 
ability, collective motion and the dynamics of heavy ion colli¬ 
sions (see, e.g. IfTTO for recent summaries). Ab initio calcu¬ 
lations of PNM with well defined theoretical errors offer ad¬ 
ditional constraints on J and L lUllTlI. Both theory and 
experiment are generally in broad agreement that L falls in the 
rather loose range 30 < L < 80 MeV, although higher values 
in particular are not completely ruled out |9| Fig. 1 shows a 
selection of experimental constraints on L, together with con¬ 
straints inferred from astrophysical observation, some of which 
will be discussed in this review ca. 

Since neutron star matter contains a low fraction of pro¬ 
tons, many inner crust and global stellar properties are sensi¬ 
tive to the symmetry energy parameters. To give a classic ex¬ 
ample, the pressure of PNM at saturation density is given in 
the parabolic approximation by Ppnm(^o)=^o^/ 3. The pres¬ 
sure at the crust-core boundary and in the outer core is domi¬ 
nated by neutron pressure so a strong correlation should exist 
between the pressure in neutron stars near saturation density 
and L. Neutron star EOSs which have higher pressures at a 
particular density are often referred to as ‘stiff’; lower pressure 
EOSs are referred to as ‘soft’. Thus ‘stiff’ EOSs at saturation 
density are associated with high values of L and ‘soft’ EOSs 
with low values of L. This fact leads to a strong correlation be¬ 
tween the radii of neutron stars and the slope of the symmetry 
energy near saturation density m. 

The crust of a neutron star is divided into two layers; the 
outer crust and inner crust. The microscopic structure of the 
crust is a crustal lattice of nuclei. The nuclei become increas¬ 
ingly more massive and neutron-rich with depth, immersed in 
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Fig. 2. (Color online) Schematic depiction of the structure of the neu¬ 
tron star inner crust from the outer crust-inner crust boundary (left- 
hand side) to the inner crust-core boundary (right hand side). Taken 
from (32. 


a gas of relativistic electrons. The distinction between the two 
layers comes from the existence of an additional component 
of interstitial fluid neutrons in the inner crust. The absence of 
free neutrons in the outer crust allow a more direct connec¬ 
tion between theoretical models and measurable exotic nuclear 
masses, and the composition and properties of the outer crust 
are substantially less uncertain than those of the inner crust 
(although uncertainty does still remain, and is connected in 
part with uncertainty in the nuclear symmetry energy ca). A 
schematic picture of the structure of the inner crust is shown in 
Fig. 2. The free neutrons in the inner crust are expected to be in 
a superfluid state, which finds some observational support from 
studying the cooling of neutron stars in X-ray transients, after 
accretion induced heating of the crust mm. At the base of 
the inner crust, competition between the inter-nuclear Coulomb 
energy and the surface energy of nuclei make it energetically 
favorable for the nuclei to form cylindrical, slab or cylindri¬ 
cal/spherical bubble shapes, a set of phases which go under the 
name “nuclear pasta” mm. Analogous changes in the mi¬ 
croscopic structure of certain terrestrial soft condensed matter 
systems such as surfactants are known, and we might expect 
the nuclear pasta phases to exhibit similarly rich and complex 
mechanical and transport properties m Because of this, it is 
useful to refer to a possible layer of the crust containing nuclear 
pasta phases as the mantle, to distinguish it from the expected 
crystalline region of the inner crust. 

The dependence on the symmetry energy of the crust-core 
transition density and pressure, and hence the crust thickness 
and extent of the nuclear pasta phases, has been well stud¬ 
ied fr8irf9ii2nii2TTi22ii^ . The transition pressure, and hence 
thickness, mass and moment of inertia of crust, correlates in 
a non-trivial way with L, i^sym, and to a lesser extent higher- 
order density expansion parameters of S{n) 11211122L The crust 
composition is also sensitive to the symmetry energy 1^ . The 
sensitivity to the symmetry energy of other crustal properties 
important for the modeling of astrophysical phenomena, for 
example the shear modulus, neutron superfiuid gap, thermal 
and electrical conductivities and entrainment of neutrons by the 
crustal lattice have yet to be examined explicitly. Nevertheless, 
we can start to examine the implications for the slope of the 


Fig. 3. (Color online). Predictions for the energy per neutron of 
pure neutron matter versus neutron baryon density obtained from our 
L=30,70 and 110 MeV EOSs using the Skyrme model SklUFSU 
mm. These are compared with calculations of Fermi gases in the 
unitary limit ||4l (SP), chiral effective field theory fS) (HS), quan¬ 
tum Monte Carlo calculations using chiral forces at leading order 
ISl (LO), Auxiliary Field Diffusion Monte Carlo using realistic two- 
nucleon interactions plus phenomenological three-nucleon interac¬ 
tions AV8-I-UIX tbllTl . and the APR EOS ISO). Figure adapted from 

Ea. 


symmetry energy L of a number of astrophysical observations 
of phenomena in which the crust is expected to play a role. We 
shall review six such observations in this paper. 

Some of the studies we review were carried out by our 
group; let us therefore give a brief outline of the EOSs we use 
in those studies. We consistently calculate the crust and core 
EOSs, the crust composition, and the transition densities and 
pressures between the inner crust and mantle, and crust and 
core, using the same model for nuclear matter. We use Skyrme 
and Relativistic Mean-Field (RMF) nuclear matter models. Two 
parameters of the Skyrme and RMF models affect only isovec¬ 
tor properties of matter and can thus systematically adjusted 
to vary the symmetry energy J and its slope L at saturation 
density while leaving symmetric nuclear matter properties un¬ 
changed II241I2511 . We construct single parameter families of 
neutron star EOSs, parameterized by L 1261127112311^ . When 
adjusting L to obtain a new EOS, we readjust J so that the pre¬ 
dicted pure neutron matter (PNM) EOS lies within the robust 
predictions of ah initio calculations at sub-saturation densities; 
thus the families of EOS contain implicitly a positive correla¬ 
tion between J and L (25 \ . Typically, we vary L in a conserva¬ 
tive range between 25 MeV and 115 MeV. We give examples 
of the predictions for PNM given by our EOSs for L = 30, 70 
and 110 MeV together with comparisons with ab initio calcu¬ 
lations, in Fig. 3. For softer symmetry energies at high densi¬ 
ties, the resulting EOS is matched onto two successive poly¬ 
tropes as described in 12711291 in order to match the constraint 
on the maximum mass of NSs of M > 2Mq OOEIl. Crust 
models are calculated using a compressible liquid drop model 
(CLDM) 1^ . The crust and core EOSs used in these works 
are available from our website http : //williamnewton . 
wordpress . com/ns-eos/ In Eig. 2 we indicate the ranges 
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Fig. 4. (Color online) Period - period derivative diagram of observed radio and X-ray pulsars. Pulsars observed in the radio band are represented 
by black dots. Anomalous X-ray pulsars are shown as the red triangles. Isolated neutron stars with thermal X-ray emission are represented by 
green diamonds, and pulsars known to be members of a binary system are shown as blue squares. Lines of constant characteristic age and 
magnetic field, inferred from the period and its derivative in a simple dipole spin-down model, are shown by the dashed and days-dotted lines 
respectively. The six observational probes of the neutron star crust that we discuss in this article are associated with the regions of diagram 
labelled (2)-(7) in reference to the section in which they are described. Data from m. 


of thickness for the inner crust and pasta phases derived from 
varying the neutron star mass between 1 and 2 Mq and the 
symmetry energy slope within the range stated above 13^ . 

Within this model, given the choice of fitting our PNM EOS 
to robust ab initio calculations, the relative thickness, mass and 
moment of inertia of the crust to the stellar total increases with 
L, whereas the relative thickness and mass of the nuclear pasta 
layers to the crust total decreases with L. 

It is useful to place the various astrophysical probes of the 
crust in the overall context of a neutron star’s life cycle. In 
Fig. 4 we plot the known periods P and period derivatives P 
of neutron stars. Standard neutron stars are believed to be born 
rapidly rotating and with magnetic fields of order 10subse¬ 
quently slowing down due to magnetic braking. Characteristic 
ages and magnetic fields, derived assuming the model of spin- 
down due to a dipole magnetic field are shown as the dashed 
and dash-dot lines respectively. Young radio pulsars appear to¬ 
wards the top left of the diagram, and over a timescale of 10^ - 
10^ years, migrate generally downwards and to the right - i.e. 
towards longer periods and smaller period derivatives (weaker 
magnetic fields as the field decays over time); eventually the 
pulsar will no longer be visible. A class of bright X-ray pul¬ 
sars, so-called anomalous X-ray pulsars (AXPs), are believed 
to be powered by extremely strong magnetic fields; such neu¬ 
tron stars are thus referred to as magnetars. They appear in the 


top right corner of the diagram (the red triangles). Such stars 
can also undergo regular bursts in X-ray and soft gamma-ray, 
as a result of magnetic field rearrangement and crust crack¬ 
ing; these objects are referred to a soft gamma-ray repeaters 
(SGRs). If a member of a binary system, in later life the neu¬ 
tron star can accrete matter off its companion star. The resulting 
transfer of angular momentum spins the star up, often to mil¬ 
lisecond periods, and the star may become visible again in the 
bottom left comer of the diagram, its emission now powered 
by accretion and thermonuclear bursts on its surface. 

A neutron star can therefore experience a rather rich life 
and exhibit a wide variety of observable phenomena; see 1351 
|36l for more detailed reviews of the neutron star ‘zoo’. Sub¬ 
stantial recent work has attempted to unite the various obser¬ 
vational classes of neutron stars under a single evolutionary 
model |37l[38l. Such diversity of observational phenomena is 
welcome as we seek a variety of independent probes of neu¬ 
tron star physics to test our models against. The regions of the 
period-period derivative diagram associated with the observ¬ 
ables we shall be reviewing in this paper are labelled (2)-(7) 
in Fig. 4 for the section they appear in. They are as follows. 
Section (2): the cooling of the neutron star in Cassiopeia A 
(this object does not have a measured P or P, so its label is 
only indicative). Section (3): upper limit of the spin period of 
X-ray pulsars (appearing as the bar in the top-right of the dia- 
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gram). Section (4): Glitches from the Vela pulsar. Section (5): 
Quasi-periodic oscillations in the light curves of giant flares 
from SGRs. Section ( 6 ): lower limit of the spin period of mil¬ 
lisecond pulsars, depicted as the bar in the lower left corner 
of the diagram. Section (7): electromagnetic pre-cursors to soft 
gamma-ray bursts, thought to result from the merger of two 
neutron stars. In each of the sections we will discuss the in¬ 
ferred constraints on L, and make clear in a series of conclud¬ 
ing caveats the main limitations of models used to obtain those 
constraints and the interpretations of observations required to 
make them. In section 8 , we shall summarize constraints on L 
and deliver some concluding comments. 


2 The cooling of the neutron star in the 
Cassiopeia A supernova remnant 

The Cassiopeia supernova is estimated to have occurred, and 
the resulting neutron star born, in the year 1680 ± 20 |[3^ . 
In 2009, thermal X-ray emission from the Gas A neutron star 
(CANS) was fit utilizing a Carbon atmosphere model BOl . giv¬ 
ing an average surface temperature (Teff) ~ 2.1 x lO^K as 
measured at the stellar surface. We note that the Carbon atmo¬ 
sphere composition is preferred on the basis that it gives an 
emitting area of order the neutron star size. Further analysis of 
data taken from the Chandra X-ray telescope, taken over the 
previous decade, indicates a rapid decrease in surface temper¬ 
ature by ^ 4% ED. A recent detailed analysis of the Chandra 
data, including all X-ray detectors and modes, provides a more 
conservative 2-5.5% temperature decline over the same time 
interval 1421 : another recent analysis concludes, however, that 
the current data is consistent with there being no cooling from 
the CANS 1431 . Even a 2% decline would constitute a rapid 
cooling phase. A definitive measurement is difficult due to sur¬ 
rounding bright and variable supernova remnant, and further 
data from Chandra over the next decade or more is required to 
settle the issue. 

Within the minimal cooling paradigm (MCP) El , the rapid 
cooling of the CANS is interpreted as the result of enhanced 
neutrino emission from neutron Cooper pair breaking and for¬ 
mation (PBF) as the core neutrons transition to a superfluid 
phase, thus providing the first evidence for stellar superflu¬ 
idity ( 145 11461 : however, alternative explanations involving in¬ 
medium modifications to the modified Urea (MU) neutrino emis¬ 
sion process in the core, neutrino emission processes in exotic 
phases such as quark matter, and heating of the crust due to 
dissipation of small-scale magnetic fields in the crust have been 
proposed 147 114 8 11491 . Core neutrons (protons) form Cooper pairs 
in the ^P 2 (^So) channel; the critical temperatures Tc for the on¬ 
set of superfluidity are strongly density dependent, and suffer 
significant theoretical uncertainty. 

The maximum value for the ^P 2 critical temperature deter¬ 
mines the age at which the PBF cooling phase is entered. Tc 
can be tuned so that the PBF cooling curves pass through the 
observed temperature of the CANS at the age of « 335 years. 
Using the APR EOS (501, the PBF cooling epoch best fit the 
data with a critical temperature of Tc = 5 — 9 x lO^K 14611451 . 
The rate of cooling during the PBF phase is determined by the 
core temperature at the onset of the PBF cooling phase, Tpbf- 



Fig. 5. Two pairs of cooling windows, together with the cooling curves 
which pass through the average value of the measured CANS tem¬ 
perature and the corresponding value of T™^^. Each plot shows the 
window with BCPs inactive (solid lines) and active (dashed lines) for 
combinations (L(MeV), M/M©, 77 ) of (30, 1.4, -13) (a) and (60, 1.4, 
- 8 ) (b)). The largest and smallest cooling rates estimated from obser¬ 
vations are indicated by the straight lines passing through the average 
CANS temperature. Taken from (33]. 


A higher Tpbf leads to a steeper cooling rate. The inferred 
^4% temperature decline is relatively steep compared to the 
predictions of theoretical cooling trajectories, and it was found 
that it could be reproduced only if the protons throughout core 
of the star had already undergone a transition to superconduc¬ 
tivity, thus suppressing the MU process in the core and raising 
the core temperature. 

In the MCP, three additional parameters affect the cooling 
trajectories of the NSs (441 : 

(i) The EOS of nuclear matter (NM). 

(ii) The mass of light elements in the atmosphere AMpght pa¬ 
rameterized here as 77 = log AMught- The thermal spectrum of 
CANS can be fit using light element masses —13 < 77 < —8 
ED More light elements means higher thermal conductivity 
and lower core temperature for a given Teff. 

(iii) The mass of CANS, constrained by atmosphere modeling 
to ^ 1.25 - 2Mq with a most likely value of 1.65M0 (5T1 . 

In the original analysis using the PBF cooling paradigm, 
this parameter space was largely unexplored; especially, the 
EOS dependence of the results was not examined. Parameter- 
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izing the EOS by the slope of the symmetry energy L, in 1331 
we examined the sensitivity of the results to the slope of the 
symmetry energy via the core EOS. 

There are additional neutrino cooling processes that have 
been postulated to operate in the pasta layers of the neutron star 
crust; if such processes are important, they will also be sensi¬ 
tive to L via the dependence of the volume of the pasta in the 
crust on L. Due to momentum conservation from scattering of 
neutrons off bubble phases of pasta, the direct Urea (DU) neu¬ 
trino cooling process might occur in the bubble phases due to 
neutron scattering |5^ . Also, neutrino and anti-neutrino pairs 
can be emitted from nucleon-pasta scattering l5^ . We shall re¬ 
fer to these collectively as bubble cooling processes (BCPs). 
The neutrino luminosity of BCPs goes as ^ lO^^Tg 

where Tg = Tcore/IO^K. The neutrino luminosity from the 
modified Urea process, assuming it to be uninhibited by su¬ 
perconductivity, goes as ^ lO^^Tg. Therefore, at core 

temperatures of order the critical temperature for the onset of 
neutron superfluidity Tg = 0.5 — 1, these neutrino emission lu¬ 
minosities are comparable. Therefore, even if superconductiv¬ 
ity suppresses the MU process, BCPs might provide additional 
cooling which will lessen the calculated cooling rate during the 
PBE cooling phase, therefore making a match to the inferred 
cooling rate more difficult. In |[33]| we tested the effect of these 
crustal processes on the inferred range of parameters for which 
the PBE cooling scenario was consistent with observation, and 
thereby obtain constraints on L through both the global cooling 
and the crustal cooling. 

To illustrate the effect of the BCPs, in Pig. 5 we show sets 
of cooling trajectories for 1.4M0 stars with two different EOSs 
(T=30 MeV (a) and T=60 MeV (b)) and envelope compositions 
(77 = -13,-8). The solid lines show the cooling trajectories 
without the inclusion of BCPs, and the dashed lines the BCP- 
modified trajectories. The upper and lower trajectories corre¬ 
spond to values of the critical temperature of OK (no ^P 2 
neutron pairing) and = lO^K (lower trajectories); each 

pair of trajectories T™^^ = OK, lO^K, forms a cooling window 
inside which the observed temperature must fall. The upper and 
lower bounds on the cooling rate inferred from the X-ray data 
are indicated by the two straight lines, which intersect at the 
average observed surface temperature at Earth (T^^), a tem¬ 
perature which takes into account the gravitational redshift of 
the emission.The trajectories which best fit {T^) are shown 
with the corresponding best-fit value of the critical temperature 
Tc. 

One can see that with the BCPs active, the temperature at 
the onset of the PBE phase is smaller and the cooling trajec¬ 
tory shallower. The effect becomes less pronounced for higher 
value of L as in those cases, the pasta phases occupy a smaller 
volume fraction of the crust. At T = 30 MeV, both cooling tra¬ 
jectories fit the inferred cooling rate (i.e. fall between the two 
straight lines); for T = 60 MeV, neither trajectory is consis¬ 
tent with the inferred cooling rate. At intermediate values of 
T, only the trajectory without BCPs will be consistent with the 
inferred cooling rate. Varying the set of parameters T, 7 /, M, 
BCPs on/off, we find the following variations lead to a lower 
core temperature and hence a shallower cooling rate which is 
harder to match to the data: 

- Inclusion of BCPs narrows the cooling window. 


- Increasing the mass of light elements in the envelope AMught 
(increasing 77 ) leads to a lower core temperature for a given 
surface temperature . 

- As M increases, the central stellar density increases and 
the fraction of the core in which the protons are supercon¬ 
ducting decreases, making the modified Urea (MU) process 
more efficient and the star’s core cooler. 

- Increasing L increases the radius R, requiring a lower Teff 
and hence core temperature to produce the same stellar lu¬ 
minosity. 

Even the 2% temperature decline over the time period 2000- 
2009 is relatively rapid; the data thus favors a relatively high 
core temperature and thus smaller values of L (smaller radii), 
smaller stellar masses M, smaller AM^ght (smaller 77 ) and 
no BCPs. With active BCPs the best fit cooling curve in the 
PBE phase is significantly less steep than without BCPs, and 
matches only the shallowest inferred cooling rate, and then 
only for the lowest values of T. As T increases beyond 50 — 
60 MeV, depending on mass, the curves become too shallow 
even with no BCPs operating. Eor each set of parameters 77 , M, 
BCPs on/off, we can find an upper limit on the value of L con¬ 
sistent with the inferred cooling rate. 

Accepting the MCP cooling model and the accuracy of X- 
ray measurements and interpretation, we obtain the following 
constraints from fitting the average CANS temperature and the 
cooling rate simultaneously: 

- Eor BCPs inactive, T < 70 MeV (35 < T < 55 MeV) for 
1.25M0 < M < I. 8 M 0 (M = I. 6 M 0 ). 

- Eor BCPs active, T < 45 MeV (35 < T < 45 MeV) for 

1.25M0 < M < I. 8 M 0 (M = I. 6 M 0 ). 

Constraint: T < 45 MeV (with pasta cooling processes), L < 
70 MeV (without pasta cooling processes) 

Caveats:. The Carbon atmosphere model is preferred because 
it gives X-ray emitting areas of order NS radii; other atmo¬ 
sphere models would shift {T^) (thus the inferred range of T), 
but not the cooling rate. Enhanced superfiuidity in the crust is 
not ruled out and would exponentially suppress BCPs, making 
our constraints closer to those obtained with no BCPs active. 
Medium modifications to neutrino cooling processes and cool¬ 
ing processes from exotic components could substantially alter 
the conclusions reached in the minimal cooling scenario. More 
recent X-ray data from the CANS is consistent with there being 
no cooling at all. 

3 Upper limit to the observed spin periods 
of young puisars 

After their birth as rapid-spinning objects, the dynamical evo¬ 
lution of most pulsars is dominated by a spin-down epoch of 
^ 10 ^ — 10 ^ years (the young - middle ages of their lives), 
caused in large part by the torque applied by their large mag¬ 
netic fields. How long a period can a pulsar spin down to dur¬ 
ing this phase of evolution? To answer this question, we must 
look to observations of X-ray pulsars rather than radio pulsars, 
because they will be free of observational biases against long- 
period pulsars that are apparent in radio observations. 
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When one looks at the periods of young to middle aged 
X-ray pulsars, the highest periods we observe are of order 10s 
(see Fig. 1). Given the lack of selection effects in the X-ray 
band, this can be taken as a physical upper limit to the period 
of a neutron star with an age < 10^ yrs. However, many of 
the longer period objects, P ^ 1 — 10s between 10^ and 10^ 
years old, also exhibit large magnetic fields IQis _ io14g 
and above, inferred from the spin-down rates. With such fields, 
these stars should spin down to periods of 100s and above on 
timescales of 10^ yrs, still within the detectable range of X- 
ray measurements. Why do we see no such long period pulsars 
with high magnetic fields? 

This inconsistency can be resolved by taking into consid¬ 
eration the decay of the magnetic field. A decay of magnetar 
strength fields to strengths of < 10^^ G on timescales of order 
10^ yrs would inhibit the spin-down of the star sufficiently that 
its period would not exceed ^ 10s 15^ . There are two main 
mechanisms for magnetic field decay in the neutron star crust 
155115611571 . At early times. Hall drift dominates the field de¬ 
cay, replaced by Ohmic decay at ages of ^ 10^yrs which are 
the ages of interest in this particular problem. At ages beyond 
^ 10^yrs, crust temperatures are cool enough that the main 
contribution to the electrical resistivity of the crust is through 
electron scattering off of impurities. A measure of the impurity 
content of the crust is given by the parameter 

Qin.p=J2xi[Zi-{Z)]^ (5) 

where the sum is over the ionic species i present in the crust, 
Xi is the fraction of species i present, is the proton number 
of the species and (Z) is the average proton number of ions in 
the crust at that particular density. Generally, Qimp ^ 1 indi¬ 
cates a relatively pure, homogeneous lattice, while Qimp > 1 
indicates an impure, heterogenous lattice. Ref. |[54l performed 
consistent magneto-thermal evolutions for high magnetic field 
neutron stars, and showed that Ohmic decay can reduce the 
magnetic field sufficiently in crust to explain the observed X- 
ray pulsar period distribution if there is a highly resistive layer 
Qimp ~ 20 — 100 at the base of the inner crust (above den¬ 
sities of 6 X 10^^ g cm“^), consistent with the existence of a 
layer of nuclear pasta which is expected to be highly disordered 
and amorphous ll58]| . Evidence for a significant pasta layer in 
the crust sets loose constraints on the slope of the symmetry 
energy L < 80 MeV 

Although the observed period distribution of X-ray pulsars 
is consistent with a significant pasta layer, it is not the only 
explanation. The structure of the inner crust of the neutron 
star is still quite uncertain theoretically. State-of-the-art semi- 
classical molecular dynamics simulations indicate the crust is 
very pure at lower densities Qimp ^ 1 1591 . Thermodynamic 
arguments, on the other hand, point towards the inner crust be¬ 
ing quite disordered throughout its extent Q 

imp 20 , also 

consistent with the period distribution I60II61I . Recently, it was 
shown that a bcc lattice is rendered unstable because of an ef¬ 
fective attractive interaction between ions induced by the free 
neutrons l62|, and that as a result the crust could assume a more 
exotic crystal structure; the effect on the resistivity is an open 
question. 

Fortunately, there is also another observational probe of the 
crustal structure. Observations of the cooling of older neutrons 


stars whose crusts are in the process of thermal relaxation af¬ 
ter a period of heating caused by accretion induced nuclear 
burning are consistent with a relatively pure inner crust, with 
Qimp ^ 1 II14II . These observations are less sensitive to the very 
deepest layers of the crust where the pasta phases are expected; 
consistency with the period distribution of X-ray pulsars sup¬ 
ports the existence of a highly resistive pasta layer. However, 
the crusts in such neutron stars may be completely or partially 
recycled, meaning that periods of accretion from a binary com¬ 
panion may have replaced some or all of the original crust, with 
the original material having sunk into the core. Recent obser¬ 
vations of cooling of a younger neutron star IGR J17480-2446 
after accretion induced heating show that the star is hotter than 
models using recycled crusts appropriate for older neutron stars 
predict 16311641 . The suggestion is that the original, cold cat¬ 
alyzed neutron star crust is more impure than the replacement. 

Despite the theoretical uncertainty, the observations are al¬ 
ready probing the structure of the inner crust in ways that give 
strong hope for tighter constraints on the crustal structure and 
thus on the symmetry energy in the future. 

Constraint: T < 80 MeV. 

Caveats: Impurity level of the pasta phases yet to be rigorously 
modeled. Some observational and theoretical evidence that the 
inner crust is quite impure throughout, removing the need for a 
highly impure pasta layer. 

4 Glitches in the Vela pulsar 

Young pulsars < 10“^yr are frequently observed to have their 
spin-down evolution interrupted by sudden increases in their 
spin frequency z/, events known as glitches. A wide range of 
glitch sizes ^ ^ 10“^^ to 10“^ and recovery timescales and 
behaviors have been reported Ea. 

In order to probe the mechanism responsible for glitches, 
one can attempt a population synthesis of all glitches observed 
to date and attempt to develop models that account for the vari¬ 
ability in the population statistically 19^ or one can examine 
a particular pulsar from which a large number of glitches have 
been observed, and attempt to account for the details of individ¬ 
ual glitches. We recap here recent progress on the latter front 
examining the behavior of glitches from the Vela pulsar. 

21 giant glitches have been observed (giant meaning ^ > 
10“^) from the Vela pulsar, and the recurrence times appear 
quasi-periodically distributed with an inter-glitch timescale of 
« 3yrs 1921 . The timescale on which the frequency increase 
occurs (the glitch rise time) has been constrained to < 40s iTSl . 

The leading class of glitch models postulate two compo¬ 
nents of the neutron star interior that for the majority of the 
time are dynamically decoupled, with one component coupled 
strongly to the magnetosphere and hence spinning down due to 
magnetic braking (this is the component whose spin frequency 
we observe), while the other component spins at a constant rate 
and hence accumulates an excess of angular momentum rela¬ 
tive to the former component (this component is often referred 
to as angular momentum “reservoir”). Occasionally the two 
components strongly couple, and the reservoir transfers angu¬ 
lar momentum to the component coupled to the magnetosphere 
and we observe a sudden frequency increase. 
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The version of this model that has gained most traction 
posits that the angular momentum reservoir consists of the su¬ 
perfluid neutron immersing the nuclei and pasta in the crust 
16811661 . The superfluid forms a two-dimensional array of vor¬ 
tices whose area density is proportional to the spin frequency; 
for the superfluid to spin down with the rest of the star, the 
vortices must move radially outwards from the rotation axis. In 
the core, the vortices are able to do this; however, neutron su¬ 
perfluid vortices in the crust interact with the nuclei there and 
become pinned in the lattice msi, preventing further out¬ 
ward motion and thus becoming decoupled from the rest of the 
star IMll^l69ll^ . As the difference between the angular ve¬ 
locity of the two components increases, the Magnus force on 
the inner crust neutron vortices grows until, when the angular 
momentum lag exceeds a threshold value, it overcomes the ef¬ 
fective pinning force and vortices unpin as a whole, coupling 
to the charged component and spinning it up. 

As a model for pulsar glitches, there are several problems 
with the current paradigm which are actively being worked on: 
there are problems with how well the model predicts the post¬ 
glitch recovery of the spin frequency to its pre-glitch behavior, 
and the details of how the glitch is initiated. However, at a more 
fundamental level, the model should be consistent with the ob¬ 
served inter-glitch timescales and the glitch sizes. Here we fo¬ 
cus on the latter in relation to the Vela pulsar, and examine 
the fundamental question of whether the crust-initiated glitch 
model can account for glitches of the size observed in the Vela 
pulsar, which sensitively depends on the EOS and, particularly, 
the symmetry energy at saturation density. 

Let us denote the moment of inertia of the angular momen¬ 
tum reservoir AI, the moment of inertia of that part of the star 
that couples strongly to the angular momentum reservoir on 
timescales of order the glitch rise time (< 405 for Vela) by I 
(which has an upper limit of the total moment of the star by 
I tot)- Deflne the parameter G as ll90ll 


G=—-> -^A 

I - |i?i 


( 6 ) 


where A is the glitch activity parameter of the pulsar, that is the 
slope of the straight line flt to a plot of the cumulative relative 
glitch size over time ll90ll . D is the average spin frequency over 
the period of time during which the glitches occurred and i? is 
the steady-state spin down rate between glitches. Then analysis 
of the rate of spin-up of the Vela pulsar caused by glitches gives 
us an observational constraint of AI/I > 1.6% I801I90I . 

It was shown that, taking / = I tot and Z\/ as the moment 
of inertia of the whole crust as an upper limit to the moment 
of inertia of the crustal superfluid neutrons Jcsf, many realis¬ 
tic neutron star EOSs can satisfy AI/I > 1.6%, and more¬ 
over that the dependence of the crustal moment of inertia on 
the EOS potentially allows constraints to be set on the symme¬ 
try energy i97iiW)ii2n . In the initial studies, it appeared that all 
but the softest EOSs predicted maximum glitch sizes consis¬ 
tent with Vela. However, the effect of entrainment of superfluid 
neutrons by the crustal nuclei was being neglected. Bragg scat¬ 
tering by neutrons off of nuclei, analogous to the similar effect 
for electrons in metals which gives rise to the band structure 
of the single particle energy levels, couples a certain fraction 
of the neutrons to the crustal lattice, thus reducing the frac¬ 



Fig. 6. (Color online). Illustration of the “snowplough” glitch mech¬ 
anism. Neutron vortices are represented by the straight black lines 
aligned with the rotation axis i7. To be pinned effectively they must be 
fully immersed in the crustal lattice, which only occurs in the strong 
pinning region of the inner crust. The red regions in the crust and core 
represent the contributions to the moment of inertia of the star from the 
charged components of the star and the core neutrons that are strongly 
coupled to them (the fraction of which is given by Yg.) The fraction of 
free crustal superfluid neutrons - that is, neutrons not entrained by the 
crustal lattice - is depicted in blue; the strength of entrainment in the 
crust is parameterized by e. 


tion of superfluid neutrons decoupled from the crust by a factor 
of, on average, ^ 5 (that is AI ^ 0.2/csf). This would ap¬ 
pear to reduce AI // too much to explain the sizes observed 
in the Vela pulsar I751l76[[70ll . However, these studies assume 
that the crust-core coupling at the time of glitch is sufficiently 
strong that the whole star spins up with the crust: I ~ /tot- 
The crust-core coupling strength is still rather uncertain, how¬ 
ever, and estimates of the coupling timescales due to interac¬ 
tions of neutron vortices and type-I or type-II superconducting 
protons ( 167119911711187117311891 ) are consistent with the possi¬ 
bility that, on the glitch rise timescale < 40s (781, only a small 
fraction of the core neutrons are coupled to the crust and there¬ 
fore / ^ /tot- This would signiflcantly reduce the ratio AI/I 
allowing it once again to satisfy the lower bound of 1.6% even 
allowing for crustal neutron entrainment at the level suggested 
by theory. 

The picture is still not complete. A recent hydrodynami- 
cal model of glitches 1841196111001 . referred to by its authors as 
the “snowplough” model, shows that vortices are pinned in the 
crust only in the region in which they are totally immersed, a 
“napkin ring” shaped annulus of in the inner crust which con¬ 
tains only ^ 10% the mass of the whole crust. This is because 
vortices that thread the core, and are only pinned at in the crust 
at their ends, can creep outwards at a rate which tracks the 
overall spin evolution of the crust, with only a very small lag 
in angular velocity. This implies that the moment of inertia of 


































W. G. Newton et ai.: Constraints on the symmetry energy from observational probes of the neutron star crust 


9 


the crust neutrons available to participate in the glitch event is 
further reduced by a factor of about 10 AI O.IAL In the 
snowplough model, / ^ /tot as outlined above, and so the ratio 
AI/I is still able to account for the observed Vela glitch activ¬ 
ity for selected EoSs 1 1001 . but the effect of crustal entrainment 
is not included. 

Recently, we performed estimates of the efficacy of glitch 
models to reproduce the observed glitch sizes in the Vela pul¬ 
sar including all the effects above that modify the critical ratio 
AI/I which we outlined above. To summarize, (i) the reduc¬ 
tion of AI by a factor of ^ 5 due to entrainment of crustal su- 
perfiuid neutrons, (ii) the reduction of AI by a factor of ^ 10 
due to the localization of the strong pinning region to an annu¬ 
lus of the inner crust, and (iii) the reduction of the moment of 
inertia of the core that is coupled to the crust at the time of the 
glitch / due to the fact that only a fraction of core neutrons Yg 
are coupled strongly to the charged component of the star then 
|[85]| . We explored the EOS and crust-composition dependence 
by using a family of SklUESU Skyrme crust and core EOSs 
(251, including the consistently calculated neutron fraction in 
the crust and core. We now summarize the specific model and 
results. 

This model is illustrated in Eig. 6. The component of the 
star that effectively acts as the angular momentum reservoir for 
the glitch - the fraction of crustal superfiuid neutrons in the 
strong pinning region that are not entrained by the crustal lat¬ 
tice, 1 — e, are shown in the blue (dark colored) region and have 
a moment of inertia I^lf ^. The charged component of the star 
includes the crust lattice plus the protons in the core and some 
fraction of core superfiuid neutrons Yg that couple to the pro¬ 
tons on timescales shorter than the glitch rise time are shown 
as the red (lighter colored) region, and have a total moment of 
inertia given by Ic. The former spins up the latter in the glitch 
event. 

What constraints do we have on Yg? In regions in which 
the protons form a Type I superconductor, neutron vortices can 
become magnetized as they entrain protons; then electron scat¬ 
tering off the magnetozed neutron vortices couples them to the 
charged component on timescales of ^ 10 — 1000s for the Vela 
pulsar I651I67I . In regions where protons are Type II supercon¬ 
ductors, they form fluxtubes which can pin neutron vortices, 
again coupling them to the charged component on timescales 
of days cmi. By comparing the above estimates of coupling 
timescales with the upper limit of the glitch rise time < 40s, 
one may infer that only some fraction of the core neutron su¬ 
perfiuid will contribute to the charged component of the star at 
the time of glitch. Because that fraction is quite uncertain, we 
take Yg to be a free parameter, but note that above estimates 
indicate that it is possible to have Yg ^ 1 18411891 . 

Let us also denote the total neutron fraction of the core at 
a given radius r by Q{r). Then the moment of inertia of the 
charged component can be expressed ifTOOl 


Ic 


Jo ‘3(0(1-^gi)]x 


e-Hr) 


Lojr) (£(r)+P(r)) 

V'l - 2GM{r)/r ’ 


(7) 


Denoting the energy density of crustal superfiuid neutrons as 
£n{r) and their pressure Pn{r), and defining 

T = ^ {Snjr) + Pnjr)) 

3 ^7 v/1 - 2GM(r)/r 

we can write the moment of inertia of the free (non-entrained) 
crustal superfiuid neutrons in the strong pinning (sp) region of 
the crust as 


/■(sp) 

^esf 
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rUn 


m^{r) 
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(9) 


where m* (r) is the mesoscopic effective neutron mass at radius 
r in the crust, values of which are interpolated from the results 
of f76ll . In order to examine systematically the effect of reduc¬ 
ing the strength of entrainment in the crust, we introduced the 
parameter e such that e=0 represents no entrainment and e=l 
represents full strength entrainment: m* ^ 1 + (m* — l)e. We 

can then calculate the ratio // as a function of the stellar 
mass M, the slope of the symmetry energy L, the fraction of 
core neutrons coupled to the crust during the glitch Yg and the 
strength of entrainment in the crust e, and confront the result 
with the observational requirement that 


> 0 . 016 . ( 10 ) 

Eor each value of Yg one can identify the value of L for 
which G > 1 . 6 % is satisfied below a given mass of neutron 
star. Eig. 7 plots those values of L vs Yg; only values of Yg 
to the left of the curves satisfy the observational constraint. 
In Eig. 7 (a-c), we show the values of L and Yg that satisfy 
G > 1 . 6 % for stars with a mass of 1.4Mq and different lev¬ 
els of entrainment e =0, 0.33, 0.67 and 1.0. Without crustal 
entrainment e = 0 , fractions of core superfiuid coupled to the 
crust on glitch rise timescales up to Yg=0.75 for L = 105 MeV 
satisfy G > 1.6% for a 1 . 4 M 0 star. Including entrainment at 
the level of e = 0.33 limits Yg <0.2; at e = 0.67, Yg <0.07, 
and at e = 1.0, Yg <0.02. In the latter case of full entrain¬ 
ment, consistency with observation is reached only for very 
high values for L >100 MeV. Because of this, we also tested 
the possibility that some of the core neutrons are also involved 
in the angular momentum transfer by extending the strong pin¬ 
ning region into the core by densities of 0.03 fm“^ Eig. 7(b), 
and 0.05 fm“^ Eig. 7(c). Extending pinning into the core by at 
least 0.03fm“^ above the crust-core transition density allows 
us to find values of L for which G > 1.6% is satisfied for any 
value of Yg (L > 95 MeV for Yg = 1.0); extending pinning 
by 0.05fm“^ allows us to satisfy G > 1 . 6 % at e = 1.0 for 
L > 57 MeV for any Yg. In Eigs. 7(d)-(f) the mass dependence 
of the results is shown taking into account entrainment at full 
strength. Eig. 7(d) shows that, of Vela was a low mass neutron 
star < 1 . 2 M 0 , the observational constraints would be met for 
L > 60 MeV for Yg <C 0 and full entrainment. 

The main results of this analysis can be summarized as fol¬ 
lows. The measured value of G from the last 4 decades worth 
of pulse timing data from the Vela pulsar generally favors a 
larger value for the symmetry energy slope L (stiffer symme¬ 
try energy) because, for this set of EoSs, a larger L tends to 
produce neutron stars with larger crusts relative to their cores, 
and hence larger AI //if AI is associated with a component 
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Fig. 7. (Color online) Lower limits to the slope of the symmetry energy L which satisfy the observational limit of G > 1.6% for a given value 
of the core superfluid neutron fraction that is coupled to the crust at the time of glitch, Yg. Panels (a-c) show results for a IAMq star with 
crustal entrainment ignored e = 0 (solid colored line), set to a strength e=0.33 (dashed line), 0.67 (dotted line) and 1 (dash-dotted line) and 
for a strong pinning region that is contained only within the inner crust nmax^ricc (a), and that extends into the outer core by a density of 0.03 
fm“^ (b) and 0.05 fm“^ (c). Panels (d-f) show results for entrainment at full strength e = 1, neutron star masses of 1.0 (solid line), 1.2 (dashed 
line), 1.4 (dotted line) and l.GM© (dash-dotted line) and for a strong pinning region that is contained only within the inner crust nmax^ricc 
(d), and that extends into the outer core by a density of 0.03 fm“^ (e) and 0.05 fm“^ (f). Observationally consistent regions of parameter space 
lie to the left of a given curve. Taken from ||8^ . 


of the crust. We specifically identify the components AI and 
I broadly with the components involved in the “snowplough” 
glitch model, and Jc respectively. Neglecting the entrain¬ 
ment of crustal superfiuid neutrons by the crustal lattice, the 
constraint on G is found to be satisfied for every value of L 
between 25 and 105 MeV, with an upper limit to the fraction 
of core neutrons coupled to the crust during the glitch event 
which increases from ^ 5% at L =25 MeV up to ^ 65% at 
1/ = 105 MeV. However, with crustal entrainment at the level 
predicted by microscopic calculations, the constraint on G is 
satisfied only with L >100 MeV and Yg ^ 1. This analysis, 
combined with previous studies 17611701 . indicate that, taking 
into account the hydrodynamic interaction between the crustal 
superfiuid and lattice, the crustal glitch paradigm is incompat¬ 
ible with observation unless the symmetry energy at saturation 
density, and hence the EoS, is very stiff. One way to resolve this 
situation is if the region of strong pinning is allowed to pene¬ 
trate some way into the core; we show that if it penetrates to 
densities up to twice the crust-core transition, then consistency 
with observations is attained for values of L down to 25 MeV 
A detailed physical account of this scenario is yet forthcoming. 


Caveats: Crust-initiated glitch paradigm still incomplete (no¬ 
tably the trigger mechanism for the glitch is quite uncertain). 
Core components could be involved in the glitch process. Fur¬ 
ther progress in falsifying models likely to involve detailed 
comparison with post-glitch spin evolution of pulsars. 

5 Quasi-periodic oscillations in the light 
curves of giant flares from Soft Gamma-ray 
Repeaters 

One observational class of neutron star, soft gamma-ray re¬ 
peaters (SGRs) exhibit regular gamma-ray fiares with total en¬ 
ergies of ^ 10^^ ergs, durations of around 0.1s, and a pattern 
of occurrence statistically consistent with self-organized criti¬ 
cal phenomena such as Earthquakes lllllll which is suggestive 
of an origin in crustal magnetic and seismic activity I109II1101 
fTT2ll . Occasionally they exhibit giant fiares with total energy 
outputs of 10^^ — 10^^ ergs. One such flare led to the discov¬ 
ery of SGRs in 1979 11061 from what is now designated SGR 
0525-66; a marginal detection of a 43Hz quasi-periodic oscil¬ 
lation (QPO) oscillation in its aftermath was reported in 1983 
Ifmll . A 2004 giant fiare from SGR 1806-20 (Ml led to a 


Constraint: L > 100 MeV 
















W. G. Newton et ai.: Constraints on the symmetry energy from observational probes of the neutron star crust 


11 


significant detection of several QPOs in the X-ray tail of the 
fiare’s light curve II11411 11 61111711 11 81 at frequencies of 18, 26, 
30, 93, 150, 626 and 1837Hz. Subsequently, analysis of data 
taken from a third giant flare, from SGR 1900-1-14, which oc¬ 
curred in 1998 11081 . revealed the presence QPOs in its light 
curve 11151 at frequencies of 28, 53, 84 and 155Hz. 

SGRs are believed to be one observational realization of 
magnetars, neutron stars with very high magnetic fields, up to 
~ 10^^G, and it is the decay of these high fields that are be¬ 
lieved to power giant fiares 1109111101 . The standard paradigm 
of giant fiares model posits that the release of magnetic energy 
in giant fiares is triggered by crust fracturing after the crust is 
stressed to breaking point by the magnetic field anchored to 
it 1109111101 . The subsequent quasi-periodic modulation of the 
post-fiare light curve is interpreted as being, at least in part, 
a result of torsional crust oscillations II 11 911 1 2211 1 2411 123 11 1251 . 
This offers a tantalizing probe of crust properties such as thick¬ 
ness and global neutron star properties such as mass and radius, 
and hence the nuclear symmetry energy, through neutron star 
seismology 112501311 . 

It has been demonstrated, for example, that varying mass 
M and radius R of a. stellar model (with fixed crust model) to 
match a particular fundamental mode frequency produces an 
orthogonal correlation between M, R, than matching the first 
overtone does. Hence, simultaneous matching of fundamental 
and first overtone has the potential to tightly constrain M and 
R I11251112611127L and hence the EOS. 

Torsional modes are designated n^i where n and I are the 
number of radial and angular nodes associated with the mode. 
Conservation of angular momentum dictates that I >2. We can 
explicitly demonstrate the dependence of the mode frequencies 
on crustal and stellar properties using a simple analytic approx¬ 
imation assuming and isotropic crust Gag. The frequency of 
the fundamental modes o^i, and of the overtones 
can be written 
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where n, / are the number of radial and angular nodes the mode 
has respectively. M, i?, R^ and A are the stellar mass and ra¬ 
dius, the radius out to the crust-core boundary and the thick¬ 
ness of the crust respectively. The symmetry energy depen¬ 
dence clearly enters through the strong correlation of R (and 
hence RA and A with the slope of the symmetry energy L. 

Vs is the shear speed at the base of the neutron star crust and 
u and A are metric fields. v‘^ = fi/p where p is the mass density 
and p is the shear modulus at the base of the crust, which, as¬ 
suming the crust has the structure of a Coulomb crystal, has a 
form extracted from molecular dynamics simulations of rnsi 
fT29lfT3()l 


Ms = 0.1106 ^ - X„)V3(^e)2. (12) 

Here, A and Z are the mass and charge numbers of the nuclei 
at the base of the crust (which occurs at a baryon density rib) 
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Fig. 8. (Color online) Fundamental (ot 2 , top) and first overtone (it 2 , 
bottom) frequencies of torsional crust oscillations as a function of the 
slope of the symmetry energy at saturation density, L, for our fam¬ 
ily of EOSs fit to ab-initio pure neutron matter calculations. Circles 
and triangles represent results with the shear modulus of pasta set to 
that of an elastic solid, while squares and diamonds represent results 
with the shear modulus set to zero. Circles and squares represent the 
cases where entrainment of crustal superfluid neutrons is neglected, 
while triangles and diamonds represent the case where it is taken into 
account. The observed QPOs from SGR flares that fall in the dis¬ 
played frequency ranges are shown by the horizontal dashed lines. 
Taken from (26). 


and Xn is the density fraction of free neutrons there. Through 
the shear modulus, the torsional mode frequency spectrum also 
depends on the microscopic structure and composition of the 
crust which is also dependent on the symmetry energy 1 ^ . 

A first study of the link between symmetry energy related 
nuclear observables and torsional crust modes using the above 
forms of the frequencies concluded that, to identify the ob¬ 
served frequencies at around 30Hz with the fundamental tor¬ 
sional mode 0^2 requires a soft symmetry energy (small L at 
saturation density) CHI. We further quantified this using our 
consistent crust-core EOSs and the same forms for the frequen¬ 
cies of the torsional modes 1 ^ to examine explicitly their de¬ 
pendence on the slope of the symmetry energy L. We also ex¬ 
amined the limiting effects of the existence of the pasta phases 
by treating the shear modulus of the pasta phases in the lim- 
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its of (i) the pasta having elastic properties similar to the rest 
of the solid crust, therefore allowing the use of the shear mod¬ 
ulus /is as given above at the crust-core transition density (a 
case we refer to as “solid pasta”), and (ii) the pasta behaves 
as a fluid with no shear viscosity, which, in this model, effec¬ 
tively lowers the density at which we evaluate the shear mod¬ 
ulus and crustal frequencies from the crust-core transition den¬ 
sity to the density at which spherical nuclei make the transi¬ 
tion to the pasta phases (a case we refer to as “fluid pasta”). 
We compared the maximal effect of pasta to the effect of su¬ 
perfluid neutrons scattering off the crustal lattice by multiply¬ 
ing the frequencies by a mesoscopic effective mass term e* = 
(1 — Xn)/[1 — Xn(mn/m*)] where m* is the mesoscopic ef¬ 
fective neutron mass ina . Physically, the crustal lattice is hav¬ 
ing to drag along a certain fraction of the dripped neutrons as it 
oscillates (in fluid mechanics parlance, the lattice entrains the 
superfluid neutrons), and the frequency is reduced. We took as 
extreme values m*/mn = 1 (no entrainment) and m*/mn = 15 
(maximum entrainment) 1771 . 

We show the results we obtained in t26i for the depen- 
dence on L of the fundamental (upper panel) and first overtone 
(lower panel) frequency in Fig. 8. Generically, the frequency 
decreases as L increases because the radius R, in the denom¬ 
inator of Eq. 13, increases with L. Circles and squares give 
the frequencies neglecting entrainment for the solid pasta and 
fluid pasta respectively, illustrating that the effect of pasta is to 
reduce the frequencies by a factor of ^ 2. The triangles and di¬ 
amonds refer to the same two cases, but with entrainment taken 
into account. Inclusion of entrainment reduces the frequencies 
by about 20-30% at most, and so is a somewhat smaller effect 
than that of the pasta phases. The horizontally dashed lines in¬ 
dicate the measured frequencies from SGRs. Allowing for un¬ 
certainties in the model, which will be discussed a little later, 
we demand only that our calculated fundamental frequencies 
fall somewhere in the range of 18 - 30Hz, the set of candi¬ 
date fundamental frequencies observed. This is obtained only 
if L < 60 MeV. If, additionally, we demand that the observed 
625Hz mode is to be matched to the 1st overtone, we then re¬ 
quire that L < 60 MeV and the pasta phases should have me¬ 
chanical properties approaching that of an elastic solid. 

One might attempt an exact identiflcation of all observed 
frequencies with theoretically calculated modes, employing the 
full spectrum of fundamental and overtones. A series of papers 
have attempted just this II341113511 1361 . A similar formalism 
for calculating the torsional modes is used, and as in our study, 
a set of crust compositions and EOSs derived as a function of 
L are employed. However, the core radius and mass is flxed in 
these studies, preventing the additional correlations of L with 
bulk stellar properties to be folded into the resulting relation¬ 
ship between L and the mode frequencies. In I134L the entrain¬ 
ment effect is neglected and the shear modulus of pasta is set 
to zero, as in our “fluid pasta” case. A similar dependence of L 
on the frequencies of the fundamental and first overtone modes 
is found, and by taking the lowest observed frequency, 18Hz 
from SGR 1806-20, to be the fundamental, they obtain T > 50 
MeV. In 11351 . using the same model but with superfluid en¬ 
trainment effects are included and the full range of observed 
frequencies from both SGRl 806-20 and 1900-^14 is simulta¬ 
neously matched to various fundamental frequencies, obtain¬ 


ing then a constraint 100 < L < 130 MeV, notably at odds 
with most experimental constraints on L. Mindful of this, CSl 
note that by excluding the second-lowest frequency observed in 
SGR 1806-20 from the analysis, one then obtains 58 < L < 85 
MeV. 

The models used in obtaining the above constraints neglect 
the effect of the stellar magnetic field, which, considering we 
believe we are dealing with magnetars, should be cause for con¬ 
cern. There are three main areas where the magnetic field can 
play an important role in the interpretation of the frequencies. 
Eirstly, the lowest observed frequencies can also be explained 
as fundamental modes or first overtones of the Alfven spectrum 
I1371I138I11391I140I (and then can potentially place constraints 
on the magnetic field geometry I141I1148I ). Secondly, the crust 
modes couple to the Alfven modes, potentially modifying the 
crust modes 1 14211143111441 or rapidly damping them out 11461 
114711148111491 . Thirdly, the magnetic field can affect the com¬ 
position of the crust itself CSol. The coupling of the crust to 
core means that the exact superfluid and superconducting states 
of matter there are highly relevant, and lead to a further enrich¬ 
ing of the possible mode spectrum and decay modes 1149111511 . 
Nevertheless, the imprint of the stellar crust and the radius and 
moment of inertia of the star on the mode spectrum will still be 
retained, and their dependence on the symmetry energy and its 
slope L cannot be ignored if we want future precision astero- 
seismology to come to fruition. 

Constraints (caveats): (i) L < 60 MeV (Theoretical frequen¬ 
cies fall in range of observed QPO frequencies; consistent crust 
& core EOSs; limiting pasta and superfluid effects crudely taken 
into account.) (ii) L > 50 MeV (Exact matching of lowest ob¬ 
served frequency to fundamental mode; inconsistent modeling 
of crust, core; fluid pasta; no superfluid effects.) (iii) 100 < 
L < 130 MeV (Exact matching of all observed frequencies; 
inconsistent modeling of crust, core; fluid pasta; superfluid ef¬ 
fects included.) (iv) 58 < L < 85 MeV (Exact matching of 
all observed frequencies except second lowest in SGRl 806-20; 
inconsistent modeling of crust, core; fluid pasta; superfluid ef¬ 
fects included.) 

General caveats: Coupling to Alfven modes ignored. Low fre¬ 
quency modes can be explained by pure Alfven modes. 

6 Lower limit to the observed spin periods 
of old, recycled pulsars 

The population of old neutron stars, evolved (or still evolving) 
through accretion-induced spin-up caused by interaction with 
a binary companion, is observed to rotate with periods up to 
milliseconds.The maximum rotation frequency of any known 
pulsar is 716 Hz 11521 . Rapidly rotating neutron star models 
constructed with reasonable EoSs tend to give theoretical max¬ 
imum spin rates of ^ 2000 Hz and above, beyond which ma¬ 
terial will be ejected from the equator of the star. This sug¬ 
gests that neutron stars are not continually spun-up by accretion 
torques until they reach their break-up frequency, but rather a 
counter-torque is encountered at substantially frequencies that 
prevents any further spin-up I153II154I . 

An origin for the counter-torque that has been intensively 
studied over the past decade and a half involves gravitational 
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radiation reaction from a class of inertial oscillation modes 
called r-modes that might be driven unstable, r-modes are anal¬ 
ogous to Rossby waves in Earth’s atmospheres and oceans as 
they have the Coriolis force as their restoring force. Below 
a certain spin period, r-modes can become unstable through 
the Chandrasekhar-Friedman-Schutz (CFS) mechanism 11551 
11561 in which gravitational radiation from the modes drives 
the oscillations to higher amplitudes. The mode becomes un¬ 
stable if the timescale for the CFS mechanism to drive the r- 
modes unstable is shorter than the timescale for the modes to 
be damped by the internal viscosity of the star; in order to ex¬ 
plain the observed cutoff, this should occur at frequencies just 
a little higher than the highest observed frequency I1571I158I 
[T591 . There are various sources of viscosity that are, in gen¬ 
eral, temperature dependent; for a given temperature, equat¬ 
ing the gravitational radiation and viscosity timescales tqr = 
Tv gives the frequency above which stars become unstable. In 
frequency-temperature space, therefore, there is a division be¬ 
tween a region where neutron stars are stable and a region 
in which they are not (the instability window). Estimates for 
the internal temperature of millisecond pulsars can be obtained 
from X-ray and Ultra-violet observations assuming accretion 
heating is balanced by core neutrino emission. If millisecond 
pulsars for which estimates of their internal temperature ex¬ 
ist are observed to fall within the instability window, then our 
physical picture is incomplete. 

The r-mode amplitude grows towards the surface of the star, 
and hence the crust is expected to play an important role in de¬ 
termining the stability of stars. This is where the role of the 
symmetry energy in determining crustal properties comes into 
play, and two studies have utilized this to set constraints on L. 
Both studies use a simple model in which the crest is assumed 
to be perfectly rigid, and at the temperatures estimated for the 
interiors of millisecond pulsars ^ 1 — 5 x 10^K, the domi¬ 
nant source of viscosity is the shear viscosity resulting from 
electron-electron scattering. 

Fet us briefly discuss the results of our study. The timescales 
(the latter in the case of shear viscosity from e-e scat¬ 
tering) are given by I1601I1611I162I . 


1 

tgr 


327rGi?2Z+2 

[( 2 / + 1)!!]2 



2/+2 pR 


f pr‘^'-+‘^dr, 
Jo 

(13) 


1 2 '+ 3 / 2 (z + i )! j2QRlp^ P f r dr 

^ 1(21 + l)\\Ii Y 7/e Jo 7c\Rc) lie 

(14) 

Here, the viscosity is evaluated at the crust-core boundary where 
the shear damping rate is greatest 11631116411161111651 . pc is the 
crust-core transition density and Rc the stellar radius at that 
density and Uc the viscosity there; it is here that the symmetry 
energy dependence on the crustal thickness enters the analysis. 
We consider the case of / = 2, with X 2 = 0.80411. The shear 
viscosity resulting from the electron-electron scattering process 
77ee = 6-0 X 10^p^T“^ (g cm“^ s“^) 1 166111671 . By employ¬ 
ing our consistent set of crust and core EoSs, we can calculate 
the frequency of onset of the CFS instability as a function of 



1 2 3 4 5 6 

T(10" K) 


Fig. 9. Frequency above which the electron-electron viscosity is in¬ 
sufficient to damp the gravitational radiation-driven r-mode instabil¬ 
ity versus core neutron star temperature for a 1.4M0 star (left) and 
2 .OM 0 (right). Results are shown for the two bounding EoSs in our 
PNM sequence: L — 25MeV and L — 115MeV, with the region in 
between shaded in grey. The positions of 4 millisecond pulsars are 
shown with two estimates of their internal temperature; one interpre¬ 
tation of observations is that the frequency of onset of the r-mode 
instability should be higher than any observed pulsar frequency at a 
given temperature in order to be consistent with those observations. 
For comparison, the neutron star break-up frequency for L — 25,65 
and 105 MeV are shown if they fall within the frequency range of the 
plot. Taken from ITTI . 


temperature for different values of L\ the results of such calcu¬ 
lations are plotted in Fig. 9 for a 1.4M0 and 2 .OM 0 star (upper 
and lower panels respectively) (271. We will focus only on the 
high mass neutron star results since ms pulsars have been ac¬ 
creting material for a large fraction of their lives. 

The solid curve shows the frequency above which the r- 
mode becomes unstable for L=25 MeV; the dotted line for 
L=\\5 MeV. The location of four neutron stars in short recur¬ 
rence time FMXBs are shown by the star symbols; the tem¬ 
perature has been estimated in two different ways I1681I169I 
leading to upper and lower bounds. Fet us focus on the object 
4U 1608. For the range of estimated temperature, the object is 
clearly in the instability window for L = 105 MeV. It should 
therefore spin down rapidly, and we should not observe it to 
have as high a frequency as we do. Therefore, to reconcile the 
model with observations, we must exclude the L = 105 MeV 
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EoS. By this argument, and adopting a conservative require¬ 
ment that the estimated temperature interval for 4U 1608 fall 
completely below the instability window, we must require that 
L < 65 MeV. 

A similar investigation, adopting the rigid crust model, con¬ 
cludes that L > 50 MeV ifTTOl . There are two important dif¬ 
ferences between this study and ours: (i) viscous damping is 
taken to act throughout the whole of the interior rather than in 
a viscous boundary layer at between the crust and the core, and 
(ii) the neutron star core and crust are not consistently mod¬ 
eled: the radius of the star is varied independently of L, when 
in fact there is a strong correlation between the two. Neverthe¬ 
less, both studies demonstrate that the position of the r-mode 
instability window is quite sensitive to the slope of the symme¬ 
try energy L. 

Let us finally outline some of the multitude of uncertain¬ 
ties that have been evinced by r-mode modeling and by careful 
observational analysis of the millisecond pulsar systems them¬ 
selves. Firstly, it is possible that the origin of the limiting spin 
period of pulsars has a physical origin not in the r-mode insta¬ 
bility but by interactions between the pulsar and its accretion 
disc I1711I1721I154I . Secondly, the assumption of a rigid crust is 
not justified; it merely sets an upper limit on the viscous damp¬ 
ing rate (and hence an upper limit to the frequency at which 
the r-mode goes unstable). However, treating the crust realisti¬ 
cally (with a finite shear modulus), it is found that the r-modes 
can penetrate the crust and therefore the relative motion of the 
crust and the core fiuid (known as the “slippage”) is reduced, 
reducing significantly the viscous dissipation there II1731I174II . 
This significantly reduces the maximum frequency a neutron 
star can be spun up to; the magnitude of the slippage and its 
frequency dependence are still sensitive to the crust thickness 
and hence the symmetry energy. Additionally, we would ex¬ 
pect such as scenario to be sensitive to the existence of nuclear 
pasta. When modeled with the inclusion of crust-core slippage, 
or with other sources of viscosity such as mutual friction be¬ 
tween superfiuid and superconducting components of the core 
and hyperon bulk viscosity, most observed systems fall within 
the instability window I1751I176L which indicates that our un¬ 
derstanding of the r-mode instability and its dependent physics 
is incomplete; the sensitivity of the latter results to the crust 
thickness has yet to be investigated, however. One way out of 
the predicament would be that unstable r-modes saturate at am¬ 
plitudes too small to be effectively damped I1771I178I . 

Constraints (caveats): (i) L < 65 MeV (consistent crust, core 
EOSs; r-mode damping via shear viscosity in a viscous bound¬ 
ary layer at the crust-core transition) (ii) L > 50 MeV (incon¬ 
sistent crust, core EOSs; r-mode damping via shear viscosity 
throughout the core) 

General caveats: Perfectly rigid crust unrealistic; r-modes likely 
penetrate the crust, reducing the crust-core slippage and hence 
the viscous dissipation there, r-modes instability might satu¬ 
rate at sufficiently low amplitudes to allow stars to inhabit the 
instability window. Other likely sources of viscosity from su¬ 
perfiuid and exotic components ignored. Limit to frequency to 
which neutron stars can be spun-up to might have physical ori¬ 
gin elsewhere, for example in star-accretion disc interactions. 
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Fig. 10. Pre-merger time tc versus frequency of gravitational waves 
emitted by the binary system /gw- Binaries follow trajectories indi¬ 
cated by the bold dashed lines for three different chirp masses M. The 
observed pre-sGRB times for the EM flares are shown by the horizon¬ 
tal dotted lines. The predicted times for the crust shattering are shown 
by the colored bars for 6 different EOSs. The EOSs have the follow¬ 
ing values for the slope of the symmetry energy L. Gs: L=93.31 MeV, 
Rs: L=86.39 MeV, SkI6: L=59.24 MeV, SkO: L = 79.14 MeV, APR: 
L=59.63 MeV, SLy4: L = 45.94 MeV (HD. Taken from 17^ . 

7 Electromagnetic precursors to short 
gamma-ray bursts 

One endpoint of the evolution of a high-mass binary system is 
a binary neutron star. In a certain fraction of such systems, the 
first-born neutron star (the primary) undergoes a common enve¬ 
lope (CE) evolution when the secondary star’s envelope over¬ 
spills its Roche lobe 11791 . During this phase, friction shrinks 
the separation of the stars significantly. The CE phase ends with 
the ejection of envelope. If the subsequent secondary supernova 
leaves behind a neutron star, the result will be a binary neutron 
star system whose separation is small enough that the stars will 
merge due to gravitational wave-induced orbital decay within 
the Hubble time. Electromagnetic (EM) and gravitational wave 
(GW) signals from neutron star mergers can provide a wealth 
of important astrophysical information including the nuclear 
equation of state (e.g., 11801 ). Large uncertainties about the 
rate of production, and hence merger, of binary neutron stars 
remain, but the predicted rates of gravitational wave detections 
for the Advanced LIGO detector tend to be or order 1-100 yr“^ 
II1811I182II . 

Gamma-ray bursts (GRBs) are extremely energetic events 
(releasing > lO^^ergs worth of ~ MeV 7 -rays) at cosmological 
distances; short GRBs (sGRBs), lasting typically < 2s, have 
been are strong candidates to be the EM signatures of binary 
neutron star mergers 1 1831 (for a review of evidence linking 
sGRBs to NS-NS mergers, see 11841 ). 

Five sGRBs have indications of precursor gamma-ray fiares 
fiares of order 1 - 10 s prior to the main burst The mech¬ 
anism for such precursor events is currently open to interpreta¬ 
tion, but one suggestion is that the emission is associated with 
the catastrophic mechanical failure of the neutron star crusts 
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induced by dynamical tidal interactions between the neutron 
stars in the immediate run-up to merger. 

Crust failure occurs when the breaking strain Cb ~ 0.1 is 
exceeded which can occur when the tidal deformation 

reaches a certain value. The failure of the crust induced purely 
by the dynamical tidal deformation of the star is unlikely to 
cause sudden catastrophic failure, with different regions of the 
crust exceeding the breaking strain at different times 11881 . If 
the released mechanical energy is deposited as heat in the crust, 
it will not diffuse to the surface in time to leave an observ¬ 
able EM signature prior to merger; the most energetic failure 
might excite seismic crustal modes which couple to the mag¬ 
netosphere, leading to another avenue for EM emission, but in 
this scenario this would occur less than Is before crustal fail¬ 
ure. 

In the model of Gsni, the dynamical tidal interaction be¬ 
tween one neutron star in its companion resonates a particu¬ 
lar global oscillation mode, and resonantly shatters the crust. 
The candidate mode with a natural frequency commensurate 
with the orbital frequency of the binary system between 1 and 
10 seconds before merger, ^ lOOHz, is a quadrupolar crust- 
core interfacial mode, and has a strong terrestrial analogue in 
a Scholte wave - a wave which travels along an elastic solid - 
liquid interface. The natural frequency of such a mode depends 
strongly on the symmetry energy, both through the dependency 
on the stellar radius and on the crust-core transition density rice 
(and the crustal shear modulus there /ice)- 

Eig. 10 shows the pre-merger time versus the frequency of 
gravitational waves emitted by the time-varying system mass 
quadrupole, /gw, which is twice the orbital frequency 11901 . 
The light-dashed lines show the pre-sGRB time for 4 precur¬ 
sor flares. The bold dashed lines show the evolution in this pa¬ 
rameter space of binary systems with a representative range 
of chirp masses M = jiMx + The hor¬ 

izontal bars represent results of calculations of the interfacial 
mode frequencies for 6 different EOSs, and are positioned at 
the frequency of the interface mode, and their heights show 
the time before merger at which the frequency of their dy¬ 
namical tides (the orbital frequency ) will come into resonance 
with the interface modes and shatter the crust. The values of L, 
given in Eig. 10, span the range ~ 45 — 90 MeV, sufficient 
variation to change to pre-merge time by nearly two orders 
of magnitude. The pre-merger times roughly correlate with L, 
with higher L values tending to give larger pre-merger times. 
It is important to note that the crust-core transition density was 
fixed at rice — 0.065fm“^ in this work, rather than allowing 
it to additionally depend on the EOS. Therefore the L depen¬ 
dence emerges from its correlation with radius R. Larger L 
gives a large star and thicker crust, leading to lower interfacial 
crustal frequencies which will come into resonance with the 
tidal forces at earlier times. 

Even accepting the scenario, there is significant uncertainty, 
but in this analysis one can claim that the best agreement with 
the observed precursor flare times comes in the region /gw ~ 
60 — 180Hz, a region which includes the SkI6, SkO and APR 
EOSs. These EOSs span a range of L « 60 — SOMeV. 

Constraint: 60 < L < 80 MeV. 


Caveats: Crust not calculated consistently with core EOS. The 
proposed mechanism for electromagnetic radiation as a result 
of catastrophic crustal failure needs a more quantitative ground¬ 
ing. The observational interpretations of EM precursors is still 
quite tentative. 


8 Conclusions 

We have reviewed recent work modeling neutron star phenom¬ 
ena in which crust physics plays a significant role with a view 
to exploring the sensitivity of the phenomena to the density 
slope of the symmetry energy L. These various models have 
been confronted with six sets of observations, taken at various 
stages of the neutron star’s life from ^ 300 yrs after its birth to 
a few seconds before its death, ~ lO^yrs later. Table 1 summa¬ 
rizes all the constraints on L obtained from these analyses. 

We start by pointing out that these diverse phenomena present 
a variety of challenges, both in the interpretation and in the 
modeling required to understand them. There are many points 
in any one of the analyses presented in which significant source 
of uncertainties can enter. Many of the observations are subject 
to possible systematic error. The theoretical models presented 
are far from complete, often neglecting important physics such 
as magnetic fields in order that they be tractable; many of them 
should be viewed as not much more than toy models. Einally, 
the derivation of constraints on neutron star physics from the 
application of theoretical models to observations is of course 
contingent on the correct interpretation of the observations, some¬ 
thing that is still open to debate in many cases. 

The sensitivity to the symmetry energy in all these cases 
emerges mainly through the sensitivity of the crustal volume 
and mass, and that of the nuclear pasta phases, to the density 
expansion parameters of the symmetry energy, particularly the 
first-order parameter L. There is also some contribution from 
the dependence of the crustal composition on the symmetry en¬ 
ergy. Various other important microscopic quantities such as 
the pasta shear modulus, breaking strain, crustal neutron su- 
perfiuid gap and entrainment, thermal and electrical transport 
properties have been taken to be fixed. Their dependence on 
the symmetry energy needs to be quantified if we are to place 
rigorous astrophysical constraints on nuclear interactions. 

Nevertheless, no matter how simplified the models, their 
demonstrated sensitivity to the symmetry is unlikely to be sig¬ 
nificantly reduced upon the inclusion of more realistic physics. 
What these studies demonstrate is that there exist a wealth of 
observations of diverse phenomena that can help constrain the 
symmetry energy in the future. Eurthermore, as constraints on 
the symmetry energy become narrower, such modeling can of¬ 
fer consistency checks and test the validity of a particular model 
or observational interpretation. Eor example, if a strict con¬ 
straint of L <70 MeV is obtained by some future experiment, 
this places serious doubts on the current pulsar glitch paradigm, 
and indicates that the model of crustal modes or the identifica¬ 
tion of crustal modes in ESI is incomplete. Thus, as well as 
astrophysical studies placing constraints on nuclear physics, in 
the future we can expect nuclear physics constraints to place 
ever tighter constraints on the theoretical models and interpre¬ 
tation of astrophysical phenomena. We can look forward to 
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Table 1. Ranges of the slope of the symmetry energy L(MeV) inferred from modeling neutron star observables 


Observable 

L (MeV) 

Specific (general) conditions/caveats 

Reference 

Cooling rate of Cas A 

:$70 

No pasta cooling processes 

EU 

neutron star 

^45 

Pasta cooling processes active and unsuppressed by crust superfluity 
(Minimal cooling paradigm; range of L contingent on atmosphere model) 


Limiting spin period of high 
magnetic field X-ray pulsars 

^80 

Magnetic field decay from highly resistive pasta layer, not 
high resistivity of an amorphous/heterogeneous inner crust 

IMI 

Vela pulsar glitches 

>100 

Full crustal entrainment, very weak crust-core coupling. 

Glitch mechanism might involve angular momentum transfer 
from core components. 

isj 

QPOs in X-ray tails of 
giant fiares from SGRs 

:$60 

Calculated frequencies fall in range of potential observed fundamental 
frequencies; consistent crust-core EOS; 
limiting superfiuid, pasta effects included 

Ei 


^50 

Exact matching of fundamental mode with lowest observed frequency 

QPO; inconsistent crust, core models; no superfiuid effects; 

C3l 


100<L< 130 

Exact matching of all observed frequency with crust modes; 
inconsistent crust, core models; superfiuid effects included 

IT351 


58 < L < 85 

As above, but with the 2nd lowest observed frequency from SGRl 806-20 

omitted in mode indentification 

(Alfven wave coupling to crust modes ignored. 

Low frequency modes could be explained by pure Alfven modes.) 

IT361 

Limiting spin-up 
frequency of 

^65 

Consistent crust-core EOS; viscous 
dissipation at crust-core boundary 

Eli 

millisecond pulsars 

>50 

Inconsistent crust-core model; viscous 
dissipation throughout entire core 

(Crust not perfectly rigid, r-mode saturation might allow stars to spin 
-up into instability window. Superfiuid, exotic shear viscosity sources 
ignored. Alternative physical mechanisms that limit spin-up are possible.) 

ITtoI 

Observed occurrence 
times of precursor 

7 -ray flares before sGRBs 

60 < L < 80 

Inconsistent crust-core EOS. Observational interpretation 
of pre-cursor gamma ray signals tentative. 

EqI 


an ever tighter interplay between nuclear physics and astro¬ 
physics. 
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